UNIQUENESS FOR SOLUTIONS OF THE TWO-PHASE STEFAN PROBLEM 
WITH SIGNED MEASURES AS DATA 



MARIANNE K. KORTEN AND CHARLES N. MOORE 



Abstract. We consider the two-phase Stefan problem ut = Aa{u) where a(u) = u + 1 for u < — 1, 
a{u) = for — 1 < M < 1, and a{u) = u — 1 for m > 1. We show uniqueness of solutions which have 
signed measures as initial data, that is, we show that if the difference of two solutions u and v defined 
on R" X (0, T) vanishes in a weak sense as 4 — » then u = v a.e. 



1. Introduction 



In this paper we show uniqueness theorems for distributional solutions to the degenerate parabolic 
equation 



in the domain K" x (0, T), for some T > 0. Here a{u) = if — 1 < u < 1, a{u) = u — 1 for u > 1, and 
a{u) = u + 1 for u < — 1. 

Equation (jl.ip is known as the two-phase Stefan problem, which describes the flow of heat within a 
substance which can be in a liquid phase or a solid phase, and for which there is a latent heat to initiate 
phase change. This allows for a "mushy zone" , that is, a region which is between the liquid and solid 
phases. In this model, u represents the enthalpy and a{u) the temperature. 

Throughout, we will consider distributional solutions to (|l.ip . that is, we consider u E Lj^^{W^ x (0, T)) 
which satisfy 



for every ip G C°° with compact support in M" x (0,T). 

For non- negative solutions of (|1.2p . a collection of techniques is available. It is possible to identify 
the growth at infinity and to show existence and uniqueness of solutions to the Cauchy problem in the 
optimal class (of growth at infinity) for measure data. See Korten [7] and Andreucci and Korten [T] . In 
this present paper, the fact that we are working with signed solutions complicates matters and we will 
need to devise a different strategy. 

Previous work of the authors (Korten and Moore [8]) provides a solution for the Cauchy problem for 
equation (II. 2|) . 

Definition 1.1. Given a Radon measure ji with /g„ exp(— c|a;p)d|/i| < oo for some c > 0, we say that 
u which is integrahle on bounded subsets of R" x (0, T) satisfies the Cauchy problem with initial data fj, 
for the two-phase Stefan problem in the sense of conservation laws if 



(1.1) 




(1.2) 





for every ip G C°°(R" x (— oo,r)) with compact support. 
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Theorem 1.2. Suppose ^ is a Radon measure which satisfies J-^„ exp(— c|a;p)(i|/i| < oo for some 
c > 0. Then there exists a solution of the Cauchy problem with initial data /i for the two-phase Stefan 
problem in the sense of conservation laws on M" x (0,T), where T = j^. 

In 8J it is shown that a solution in the sense of conservation laws is a solution in the sense of 
distributions. Consequently, we have: 

Theorem 1.3. 8J Suppose ji is a Radon measure satisfying Jjj„ exp(— c|xp)(i|/i| < oo for some c > 0. 
Then there exists a solution of (jl.2p on M" x (0, T), T = which has initial value fi in the sense that 
linit^o /jjn u{x, t)^p(x)dx = J^^ 'il){x)dfj, for every function ^ € C°°(R") of compact support. 

Note that a solution in the sense of conservation laws is integrable on bounded subsets of R" x (0, T), 
whereas a distributional solution is required only to be integrable on compact subsets of K" x (0, T); for 
u which are integrable on bounded subsets of R" x (0, T), the two notions of solution are equivalent. 

In other previous work, the authors have discussed regularity of solutions. In Korten and Moore 
[9] the authors show intrinsic energy estimates for signed local solutions to (|1.2p . These estimates do 
not involve initial or boundary data. Subcaloric estimates are then used to show that if m G Lf^^ is a 
solution of the two-phase Stefan problem, then a{u) is locally bounded. It then follows from a theorem 
of Caffarelli and Evans 4J (or similar results by Sacks |T^, Ziemer [T3|, or DiBcnedetto [6]) that a{u) is 
continuous. 

Our purpose in this paper is to show the following uniqueness theorem. 

Theorem 1.4. Suppose u and v are solutions of (II. 2p on R" x (0,T) which belong to L°°(R" x {e,T)) 
for every e > 0, and which satisfy 

[ [ {\u{x,t)\ + \v{x,t)\)e-''^''^^ dtdx <oo, 

JR" Jo 

for some c > 0. If for every ip e C°°(R") of compact support we have 

lim / {u{x,t) — v{x,t))ip{x)dx = 0, 

then u{x,t) = v{x,t) a.e. onR"x(0,T). 

The hypotheses of the theorem essentially say that u = v initially in the sense of measures. The 
novelty of our result is that it allows for signed solutions. Pierre [TT] demonstrated a similar uniqueness 
result for nonnegative solutions u to equations of the form ut — /S.Lp(u) = 0, where is assumed to be 
nondecreasing, locally Lipschitz, and with ip{Q) = and m G Li(R" x (0, T)) n L°°(R" x (e, T)), for every 
e > 0. Bouillet [5] considers solutions u (possibly of changing sign) of Ut — Aa(u) = 0, which satisfy 
the same growth conditions as in our theorem, with Lipschitz a(u), and obtains uniqueness under the 
stronger hypothesis that \\u{-,t) — v{-,t)\\i^i (^n) ^ 0, as i ^ 0, when n > 2. For n = 1, Bouillet has 
obtained the result above. Related to the two-phase Stefan problem is the porous medium equation 
Ut = Aw™, TO > 1. This has been studied extensively by many authors, including uniqueness results for 
nonnegative solutions. 

We do not know optimal conditions under which a solution u belongs to L°°(R" x (£,T)) for every 
£ > 0. Korten and Moore \9\ show that the assumption u G L'j^^(^^ x (0,T)) imphes a{u) is continuous, 
hence u is locally bounded. For nonnegative solutions, Korten [7] shows that if u G L\^^ is a solution of 
(jl.2p then a{u) is continuous, but this is an open question for solutions which change sign. 

We remark that on a fixed compact set K C R", the functionals Tt{(p{x)) = Jj^{u{x, t)—v{x, t))(p{x)dx, 
for if supported in have limt_,o 7* ((/?) = 0. In particular, Tt{ip) is bounded for fixed ip, so by the 
principle of uniform boundedness, ||Tt|| = \\u{-,t) — w(-, is uniformly bounded. 

The paper is structured as follows: In section 2 we establish some technical lemmas which will be used 
in the proof of Theorem 1.4, which we prove in section 3. Section 4 contains some further discussion. 
Throughout, the letter C will denote a constant which may vary from line to line. 

Marianne Korten is supported by grant #0503914 from the National Science Foundation. 
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2. A FEW TECHNICAL ESTIMATES 



The estimates in this section involve elementary computations with the Gaussian kernel and the 
maximum principle. Throughout, B{R) denotes the ball centered at of radius R. Some of the ideas in 
this section are taken from the previously mentioned work of Bouillet [2] . 

Lemma 2.1. Consider [1, R] C M and on [1,R\ x [0,00) let w{x,t) he the solution to wt = Aw, with 
initial condition w{x,0) — 0, and boundary conditions = 1 and w{R,t) — for all t > 0. Then 

(i) ^>0 for all x £ (1, R),t>0 

(a) \^{R,t)\ < exp(^^) whenever <t <T if R is sufficiently large. 
Proof. Observe that ^ satisfies the heat equation with boundary values 



for 



t > and initial values ^{x,0) > 0. Then (i) follows immediately from the maximum principle. For 
(ii) consider the function w^x^t) on (— oo,i?] x (0, cxd) which satisfies Wt — Aw, w{x,0) 4 if a; < 1, 
w{x, 0) = if 1 < x < i?, and ■w{R, t) = for all t > 0. This can be given explicitly as: 

w{x, t) = 



'A-Kt J- 



e ds 



Now 



w{l,t) 



/A-Kt 



' ds 



'Ant J- 



-ds 



e-'^'du. 



If R is chosen sufficiently large, then for i < T the last integral is bounded by Thus, for t < T, 
'w{l,t) > 1. So by the maximum principle, w{x,t) < w{x,t) for 1 < x < R and < t < T. Since 
w{R,t) = w{R,t) = then \^{R,t)\ < \^{R,t)\. Explicitly differentiating we find 



ox 



/A-Kt 



s — R -(a-fi)^ 
e ft ds 



2-y-t 2y 2 

— -=e " ay 
30 vt 



A -a-R)^ 
e 



TTt 



If R sufficiently large, this last quantity is < -^e s* for all < t < T. 



□ 



Lemma 2.2. Suppose < f < I is supported in the ball B{1) C K", d{x,t) is a smooth, bounded 
function on B{R) x (0, T), and k < d{x,t) < 1 for some k > 0. Suppose R > 1 is sufficiently large (as 
large as required for (ii) in the previous lemma) and let h{x,t) satisfy 



Then for {y,t) with \y\ 
(2.1) 



dh 

d{x,t)Ah{x,t) = -Q^{x,t) 

h{x,0)^fix), 
h{x,t) = 0, 

R,teiO,T), 

dh 



on B{R) X (0,T), 
X G B{R), 
{x,t) e dB{R) X (0,T). 



dn 



{y,t) 



< C exp 



-R' 



Here n denotes the outward normal to the ball B{R) and C is an absolute constant (in particular, it's 
independent of R, k and d(x,t)). 

Proof. Consider (y, t) where y is of the form y — {R, 0, . . . , 0). By rotation, it suffices to show (12. ip for 
such {y,t). 

Set D = B{R) n {x e R" : > 1}. On _D x [0, 00) set W{xi, . . . , Xn,t) = w{xi,t), where w is as in 
Lemma Then AW = Wt. By the maximum principle, \h{x, t)| < 1 on B{R) x (0, T), so in particular, 
W{x,t) > h{x,t) whenever xi = 1. On [dB{R) n {{x,t) : xi > 1}) x (0,T), W{x,t) > = h{x,t). 
Consequently, h{x,t) - W{x,t) < on dD x (0,T). When t = 0, x e D, W{x,t) = = h{x,t). Thus, 
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h ~ W < on the parabolic boundary of D x (0, T). Furthermore, on D x (0,T), since d < 1, and 
AW > by Lemma [2lTi). then 

dA(/i - VK) = dA/i - dAVF > dAh - AW ^h- Wt- 

Therefore, by the maximum principle, h{x, t) < W{x, t) on Dx (0, T). In particular, for x — {xi, 0, . . . , 0), 
Kxi < R, h{x, t) < W{x, t). This combined with = h{R, 0, . . . , 0, i) = W{R, 0, . . . , 0, gives 



dh 
dn 



(i?,0,...,0,i) 



dh 

dxi 



(i?,0,...,0,0 



< 



dW 

dxi 



(i?,0,...,0,t) 



dw 
dx 



< Cexp 



8t 



□ 



Lemma 2.3. Let u be a solution of (II. 2p such that u E L°°(M" x (e,T)) for every e > 0. Then for a.e. 
<ti <t2 <T, and a.e. R> 0, 



/ u{x,t2)(l){x,t2)dx = / u{x,ti)(j){x,ti)dx 

J B{R) J B(R) Jt 

rt2 



dB(R) 



a(u) — dadt 
on 



f " f d(j) 
I I u— — \- a{u)A(l)dxdt 
Jti Jb{r) ot 



for every 4>{x,t) G C^{B{R) x {ti,t2)) H C{B{R) x [^1,^2]) which vanishes on dB{R) x [ti,t2]. Here n is 
the outward normal to dB{R). 

Proof. Let (p{y,s) be a smooth radial function of compact support on and for m = 1,2,3,... 

set ipmiy,s) = mip{my,ms). For {x,t) S M" x (0,r) and m sufficiently large (depending on {x,t)), 
Vm{x — y,t — s) is a test function and thus 



(2.2) 



'u{y,s)^^{x -y,t~ s)+ a{u{y, s))ALp,n(x -y,t~ s)dyds = 0. 
at 



Set Um — fm * u and Wm = frn * Q:{u) . Then on compact subsets of R" x (0, T), if m is sufficiently large, 
we can rewrite (|2.2p as the pointwise equality —-^Um + Awm = 0. 

Since u and a{u) are locally integrable, then a.e. point of M" x (0,T) is a Lebesgue point of both 
functions. By Fubini's theorem, for a.e. t, < t < T, a.e. point of M" x {<} is a Lebesgue point of u 
and a{u). Similarly, for a.e. R, a.e. point of dB{R) x (0,T) is a Lebesgue point of m and a{u). Consider 
< ti < t2 < T, with the property that a.e. point of R" x {ti} and R" x {^2} is a Lebesgue point of u 
and a{u), and consider R > with the property that a.e. point of dB{R) x (0, T) is a Lebesgue point of 
u and a{u). Then for sufficiently large m, and with </> as in the hypotheses, we have: 



= 



B(R) 



dUrn 

dt 



(j)Aw,n dxdt 



Ujn{x,t2)(j){x,t2)dx - I U„i{x,ti)(t){x,ti)dx - 

B(i?) JBiR) Jt 



B(R.) 



—dxdt 
ot 



—Wmdadt — 
on 



WrnAcj) dxdt. 



Iti JdB{R) on JB(R) 

Let TO ^ 00 and rearrange to obtain the conclusion of the lemma. 



□ 



In all subsequent uses of this lemma, we will assume that the ti and t2 as well as the R are chosen in 
the sets of full measure for which the above formula is valid. 
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3. The proof of the Theorem 

Fix 6 (a;) G C°°(M") of compact support. Consider to with to < min{^,T}. We will show that for 
a.e. such to, Jg^„{u{x,to) — v{x,toj)Q{x)dx = 0; that this is then true for a.e. to, < to < min{-^,T} 
and such implies u — v a.e. on R" x (0, min{^, T}). If < T, then note that the hypotheses will 
now hold on ( j^^, min{ T}) and repeating the argument gives that u = v a.e. on (0, min{ T}). 
Continuing, eventually we obtain u = v a.e. on R" x (0,T). 

Without loss of generality we will assume that 8 is supported in B{1), Q > 0. Using the differential 
equation, we will express J^ri{u{x,to) — v{x,to))Q{x)dx in terms of several other quantities. By appro- 
priately choosing certain parameters, we will show how to make each of these quantities as small as 
desired. Some of the choices of parameters will depend on the choices of other parameters, that is, the 
order in which they are chosen is critical. At the end of the proof, we explain the exact order in which 
to choose these parameters. 

We follow a technique of Oleinik jTO] . Define 



[" a(u{x^t))—a{v(x^t)) 



if u{x, t) 7^ v{x, t) 
if u{x, t) = v{x, t) 



Then < c{x,t) < 1 a.e. on x (0,T). 

Suppose R is much larger than 1, to be chosen precisely later. For m 
regularization of c{x, t) V ^ which satisfies 



1,2,... let Cm{x, t) be a 



lim 

"'-♦""Jo JB(R) 

For m — 1,2, . . . let (j)m be a classical solution of 

0t + Cm^4> = 

0(x,to) = 6(x), 
= 

< lie 



■ dxdt — 0. 



on B{R) X {0,to), 
X e B{R), 
{x,t) e dB{R) X (0,io]- 



By the maximum principle, ||(/'m||oo 
by Lemma 12.31 we may write 



B(B) 



{u{x,to) — v{x,to))Q{x)dx 



Suppose Q < 5 < to, where 5 is to be chosen later. Then 

{u{x, 5) — v(x, 5))4>rn{x, S)dx 



dB(R) JS 



to 



dn 



[a{u) — a{v)] dtdfj - 



B{R) JS 



A(/)„i [a{u) - a{v)] 



dt 



v\ dtdx 



= 1 + 11 + III, 

each of which we investigate separately. 

To estimate III, we use — CmAt/fm = §[(t>m to obtain 



I///I 



B{R.} JS 



A4>m [a{u) - a{v)] - CmA4>m [u - v] dtdx 



< \\u — V 



A0„C [u- v]- Cjn'A4'n 



v] dtdx 



L~(B(i?,)x(A-,to]) 



< ||u - v\\L°-{B{R}x(S,to]) 




\A4>m\\c - Cm\dtdx 



i\A4)rn\'^dtdx 




dtdx 
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Multiply the equation ^( 



(3.1) 



B{R) 



\V<j),n{x,S)\^dx + 



Cm^'pm = by A(/>m and integrate by parts to obtain 

2, 



S JB(R) 



\A(t)jn\ dxdt 



1 



BiR) 



|Ve(a;,io)| 



and thus /// ^ as m ^ oo (for R and 6 fixed). 

To estimate //, we use Lemma |2.2[ with t replaced by ~ t: 



II < 



< 



dB(R) JS 



dB{R) J 5 



ta 



exp 



dn 

( 



< exp cR 



\S{to-t) 
R^ 



{\a{u{x,t))\ + \a{v{x,t))\) dtda 

exp(c|a;p) (|a(u)| + \a{v)\) exp{-c\x'^\)dtda 

■ T 



JdB(R) Jo 



i\a{u)\ + \a{v)\) e-Kp{-c\x \)dtda 



Since {\a{u)\ + \a{v)\) exp{—c\x'^\)dtdx — M < oo, there exists a set G C [0, oo) with the property 

that G n [L, oo) has positive measure for any L > 0, such that for R £ G, 



dB(R) Jo 



(|a(u)| + \a{v)\)expi-c\x''\)dtda < 



Furthermore, exp (^cR 



a R^ 

8(to-t) 



> as i? ^ oo since < §7:, and hence, to — t < Therefore, we 

may choose R sufficiently large (in G), and independent of S and m so that // is as small as desired. 

We would now like to estimate / = Jgi^j^.^{u{x, S) — v{x, S))(f>mix, S)dx. We use a variation of the 
strategy so far. Let qm{x,t) be the solution to 

qt + Aq^O on B(i?) x (-00, (5), 

q{x,S) = (j>,nix,S), X e B{R), 

q{x,t)=0 on dB{R) X {-00, 6). 

Let < 7 < (5. Then by Lemma [ 



B{R) 



{u{x, 5) — v{x, 5))4>rn{x, S)dx 



{u{x, 7) - v{x, 7)) qmix, "f)dx 



ldB{R) J-t 
= /l + /2 + /3 



dqm 
dn 



[a{u) — a{v)]dtda + 



B(R.) J 7 



Aqm[a{u) — a{v)] + "^b^ ~ v]dtdx 



We first estimate I3 in a similar fashion to our estimation of ///. Since — —Aq, 



< 



B{R} 



Aqm [(x{u) — a{v) — (u — v)] dtdx 



< 2 



B{R) J-f 



I Aqm I dtdx 



2\ f f Idtdx] if [ \Aqrrfdtdx] = 2y'\B{R)\V5 I [ [ \Aqm\'^dtdx 

\JB{R.)J-t I \JB{R)Ji I \JB{R)J-i 



J \JB{R.} 



B(R) J-i 



Similar to (j3.ip , multiply the equation = + Aqm by Aqm and integrate by parts to obtain 



\! \Vqmixn)\^dx+ [ [ \Aqmfdxdt:^l [ \Vqmix,S)fdx 
^ J B(R) J B(R) J-f ^ J B(R) 

= \ I |V0„(x,<5)pdx < i / \VQ{xM)?dx 
^ Jb(r) ^ Jb(r.) 



where we have used (EH). Thus, I/3I < Gy/\B{R)\VS. 
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Our estimation of I2 is similar to that of //. We claim < Cexp( g^'^^^^j ) for (x, t) e dB{R) x 

(0, S). This follows from a variation of Lemma l2.2l Consider that lemma and its proof. Let hhe a solution 
of the equation as in the statement of the lemma; then as concluded there, i)| < C exp(^-) for 

= R and < i < T. Suppose however, < Ti < T is fixed, and define r{x,t) as the solution of 
Ar = rt on B{R) x {Ti,T), r{x,t) = on dB(R) x [Ti,T), and r{x,Ti) = h{x,Ti), for x G B{R). 
Then, following the proof of Lemma [2.21 (with the notation there) we have that h{x,t) < W{x,t) on 
D X (0, r), so in particular, r{x, Ti) = h{x, Ti) < W{x, Ti) for x £ D. Furthermore, reasoning exactly as 
before, r{x, t) < W{x,t) on dD x (Ti, T). So by the maximum principle, r{x, t) < W{x,t) on Dx (Ti, T). 
Continuing as in the proof of Lemma l2.2l we conclude | ^lii^OiiiiiLl) | < (j gxp(^_)^ and hence by rotation, 
l^^^i^l < Cexp(^) for all {x,t) G dB{R) x {Ti,T). Applying this to qm{x,to - t) yields the desired 
estimate |^2:|l£2£)| < C'exp(g^=^) for G dB{R) x (0,(5). 

With this estimate in hand, the estimation of I2 follows exactly the same steps as the estimation of 
//. Consequently, we may choose R sufficiently large, independent of 5, 7 and to, so that I/2I is as small 
as desired. 

We finally estimate Ii. 

v{x,j))qmix,-f)dx = 

v{x,j)){q,n{x,-f) ~ qjn{x,0))dx + / {u{x,'j) - v{x,j))qm{x,0)dx. 

JB(R) 

If 7 is small, the next to last integral is small since ||m(-, t) — «(•, t)\\L'^(B{R)) ^i'^ bounded (see the remarks 
after the statement of the theorem) and the fact that qm{x, t) q-mix, 0) uniformly as t ^ 0. The last 
integral is small if 7 is small by hypothesis. 

Let us explain the order in which various constants are chosen. First R should be chosen so that // 
and I2 are small. Then 6 should be chosen so that is small. Then to should be chosen so that /// is 
small. Then 7 should be chosen so that Ii is small. This forces Jg^^j (u(a;, to) ~ v{x,tQ))Q{x)dx to be 
as small as desired, which completes the proof. 




4. Further remarks 

With only slight modification, the proof of Theorem II .41 can yield a more general statement. 

Theorem 4.1. Suppose a ; R — > M is nondecreasing and Lipschitz, and that there exists a number a > 
so that a{u) — au is bounded. Suppose u and v are solutions of (|1.2p on M" x (0,T) which belong to 
i°°(M" X (e,T)) for every e > 0, and which satisfy {\u{x,t)\ + |z)(a;, t)|)e~'^l^l dtdx < 00, for some 

c > 0. If for every ip G C°°(K") of compact support we have liming J^n(u{x^t) — v{x,t))ip{x)dx = 0, 
then u{x,t) — v{x,t) a.e. onM"x(0,r). 

Neither Theorems 11.41 or 14.11 allow for a which arc only locally Lipschitz such as a{u) — sgn{u)\u\"^ 
in the porous medium equation. Benilan, Crandall and Pierre [3] have shown uniqueness, in the sense 
of distributions, for solutions u G C{[0,T) : Lj^^{M.^)) of the porous medium equation which satisfy a 
certain growth condition. But for the porous medium equation, uniqueness for signed solutions with 
the initial data taken as measures, remains an open problem. See Daskalopolous and Kenig [S] for a 
discussion of uniqueness results for the porous medium equation and this open question. 

An interesting first step toward more general uniqueness results would be to show the main theorem in 
the case when a is only assumed Lipschitz (or even locally Lipschitz) and nondecreasing. (With possibly 
different growth conditions on the solutions.) These are the assumptions in the theorem of Pierre |llj on 
uniqueness for nonnegative solutions, and the theorem of Bouillet [2] on uniqueness for signed solutions 
which assumes convergence as t ^ 0. 
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